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Abstract 

We illustrate the use of computational methods for the study of baryon 
non-conservation in high energy electroweak processes. 


1 Introduction 

Since the advent of the modern scientihc method, knowledge of the physical 
world has advanced through the interplay of experimental observation and theoreti¬ 
cal analysis. Correspondingly, the terms experimental physics and theoretical physics 
have been used to characterize these two basic methodologies of scientihc research. 



The introduction of computers, however, has lead to the establishment of a third 
fundamental methodology and the term computational physics has been coined to 
denote all those investigations where computers are used in an intrinsic manner to 
unravel properties of complex systems. The methods of computational physics do 
not offer a replacement for either experimental investigation or theoretical analysis. 
Indeed, although computer modeling shares many of features of actual experiment, 
clearly the study of computer simulated phenomena can never replace the observation 
of the actual world. In a more subtle way, computational physics cannot displace the¬ 
oretical physics either. The insights into the phenomena provided by the theoretical 
analysis constitute the necessary platform over which a computational investigation 
can be launched and are ultimately required to interpret the latter’s results. Thus 
theoretical physics and computational physics complement each other and, at their 
best, proceed hand in hand. 

Because of the large computer power required by some numerical investigations, 
it is frequently thought that computational physics is accessible only to a limited con¬ 
stituency of researchers fortunate enough to have access to powerful supercomputers. 
This belief is wrong, on two accounts. First, not all computational investigations re¬ 
quire very large resources. In many instances, ingenuity more than powerful hardware 
is the key to a successful utilization of the computer. Second, the amazingly rapid 
pace of technological development in the held of computers is making these tools ever 
more accessible on a world wide scale. For these reasons, the International Center for 
Theoretical Physics has been offering colleges on computational physics on a rather 
regular basis since the hrst conference and college it held in 1986. 

In particle physics some of the most successful computational applications have 
taken place in the area of lattice gauge theories, where computer simulations of quan¬ 
tum huctuations have been used to evaluate several non-perturbative quantities of 
primary importance. In particular, computational investigations of QCD have pro¬ 
vided evidence of quark conhnement as well as the means of calculating hadronic 
observables such as meson and baryon masses, some weak matrix elements, the de- 
conhning temperature and the value of as- However, the span of computational 
particle physics investigations is by no means limited to lattice QCD, and in these 
lectures we will illustrate another challenging application of numerical techniques to 
particle phenomena where computational and theoretical methods complement each 
other. We will see, indeed, that the theoretical analysis of the phenomenon provides 
the basis for the computational investigation, which in turn can produce information 
that goes well beyond the reach of purely analytic techniques. 

The problem we will study is the possible occurrence of baryon number viola¬ 
tion in high energy electroweak processes. In perturbation theory, baryon number 
is strictly conserved. But, as has been well known since the pioneering work of ’t 
Hooft[^, the axial vector anomaly implies that baryon number is not conserved in 
processes which change the topology of the gauge helds. Baryon number violating 
amplitudes are non-perturbative, and viable methods of calculation are scarce. Basi- 
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cally there are two ways of getting at non-perturbative information in quantum field 
theory. One can use either semi-classical techniques or direct lattice simulations of the 
quantum fluctuations. Unfortunately, theories with small coupling constants are not 
suited for the latter, so the electroweak sector of the standard model lies beyond the 
reach of direct lattice calculations. This means that semiclassical methods presently 
offer the only way to study baryon number violating electroweak processes. 

Electroweak baryon number violation is associated with topology change of the 
gauge helds. Classically, gauge held conhgurations with different topology (i.e differ¬ 
ing by a topologically non-trivial gauge transformation) are separated by an energy 
barrier. The (unstable) solution of the classical equations of motion which lies at 
the top of the energy barrier is called the sphaleron^. At energies lower than the 
sphaleron energy, topology changing transitions, and hence baryon number violation, 
can only occur via quantum mechanical tunneling. Under certain circumstances, 
semiclassical methods can be used to approximate these tunneling rates. The rele¬ 
vant solution of the Euclidean equations of motion which describe such tunneling is 
known as the instantonPJ. 

All transitions that change topology involve helds of order 1/g and actions 
S = So/g"^, where So is a coupling constant independent action reexpressed in terms 
of rescaled helds gA^{x). As a consequence, all tunneling amplitudes contain a barrier 
penetration factor exp(—S' 0 / 5 '^): 

A ~ exp(-^o/£/^), (1) 

where So is a numerical factor of order one. The appearance of the square of the elec¬ 
troweak coupling constant g in the denominator of the exponent in Eq. |l] has three 
important implications. First it tells us that the phenomenon is non-perturbative, in 
the sense that A has an essential singularity at 5 ^ = 0 and thus cannot be expressed in 
terms of a perturbative expansion in powers of g. Second, the fact that the actual nu¬ 
merical value of g is very small indicates that non-perturbative semiclassical methods, 
based on saddle point expansions of path integrals around solutions of the classical 
equations of motion, are likely to produce reliable results (to better understand this 
point, remember that in units in which h is not set equal to 1, Eq. |l| takes the form 
A ~ exp[—S'o/(h 5 f^)] so that small g and small h are equivalent). Third, the small 
value of g also makes these processes apparently irrelevant because the associated 
transition rates, proportional to |Ap, turn out to be abysmally small. 

This state of affairs changed a few years ago when Ringwald ^ and later 
Espinosa [^] noticed that a summation of the semiclassical amplitudes over hnal states 
gives rise to factors which increase very rapidly with increasing energy. This might 
lead to a compensation of the suppression factor in Eq. for energies approaching the 
energy of the barrier, i.e. the sphaleron energy Egph- Intuitively, one might expect 
the tunneling suppression factor to become much less severe as the energy approaches 
the energy of the barrier. In particular one might expect it to disappear altogether 
for E > Egph, i.e. in the region where the topology changing processes are classically 
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allowed. Investigations have indeed confirmed that this is precisely what happens 
in high temperature electroweak processes 0 . As the temperature approaches Egph 
the barrier-penetration suppression factor becomes progressively less pronounced, and 
electroweak baryon number violation becomes unsuppressed altogether for tempera¬ 
tures comparable to the sphaleron energy. The situation is, however, much less clear 
for high energy collisions. Phase space considerations are more subtle, and simply 
because one has enough energy to pass over the barrier does not guarantee that one 
does. The problem is that in high energy collisions the initial state is an exclusive 
two particle state, which is difficult to incorporate in a semiclassical treatment of the 
transition amplitude. 

A possible remedy to this situation has recently been proposed by Rubakov, Son 
and Tinyakov0, who suggested that one considers inclusive initial coherent states, 
but constrained so that energy and particle number take hxed average values 


E = 


( 2 ) 


^=- 2 - 


( 3 ) 


In the limit g ^ 0, with e and u held hxed, the path integrals giving the transition 
amplitudes are then dominated by a saddle point conhguration which solves the clas¬ 
sical equations of motion. This permits a semiclassical calculation of the transition 
rates. Information on the high energy collision processes can then be obtained from 
the limit i/ —0. While this limit does not strictly reproduce the exclusive two-particle 
initial state, under some reasonable assumptions of continuity it can be argued that 
the corresponding transition rates will be equally suppressed or unsuppressed. 

In these lectures we will not reproduce the derivation of the saddle point equa¬ 
tions, which would form the topic of an extended set of lectures in itself. (Indeed, 
Rubakov presented such lectures at the 1992 ICTP Summer School in High Energy 
Physics and Cosmology, and the conversations that one of us held with him then 
stimulated the investigation we describe here.) Rather, we will start from these equa¬ 
tions, which we will of course recapitulate, and describe the computational techniques 
used to solve them and the progress we have made in this direction. For the actual 
derivation of the equations the reader should consult Ref. 0. 

In the next section we illustrate the general properties of topology changing 
evolution of the classical helds. For simplicity we hrst consider the 2-dimensional 
Abelian Higgs model. Then we examine the 4-dimensional SU{2) Higgs model, but 
restricted to the spherical Ansatz to obtain a computationally tractable system. In 
Section 3 we investigate the properties of topology changing processes above the 
sphaleron barrier, i.e. in the classically allowed energy domain (see also Ref. |^). 
And in Section 4 we hnally describe the equations introduced by Rubakov, Son and 
Tinyakov]^] and the computational methods required to solve them. 
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2 Topology changing field evolution 

The 1+1 dimensional Abelian Higgs system is defined in terms of a complex 
matter field ^(a;) and an Abelian gauge potential A^{x) with action 

S^Jdx^ + (D,4,)‘D‘‘4,-MW - 1)"} . (4) 

where the indices run over 0 and 1, = d^A^ — dyA^, D^cj) = d^<p — tA^<p and many 

inessential constants have been eliminated by a suitable choice of units. 

The most important feature of this system is that the vacuum, i.e. the configu¬ 
ration of minimum energy, occurs for non-vanishing 0, indeed, with our special choice 
of units for |0| = 1. Since this does not specify the phase of 0, there is not a unique 
vacuum state, but rather multiple vacua. Still, because of gauge invariance one must 
be careful in regard to the physical significance of the phase of 0. A local variation of 
the phase of 0 can always be undone by a suitable gauge transformation. And since 
gauge equivalent configurations must be considered physically indistinguishable, local 
variations of the phase of the matter field do not lead to different vacua. However, 
variations of the phase of 0 by multiples of 27i (as the coordinate spans the en¬ 
tire spatial axis) cannot be undone by a local gauge transformation, and thus define 
topologically distinct vacuum states. These vacua differ by the global topological 
properties of the field configuration. The condition |0| = 1 restricts the values of the 
matter field to the unit circle (in the complex plane). If we demand that 0 takes fixed 
identical values as x^ —> ±cx) (a condition we later relax), then the number of times 0 
winds around the unit circle as x^ spans the entire real axis is a topological invariant 
(the winding number) and characterizes different topologically inequivalent vacuum 
states. 

Figures la-c illustrate three contours traced in the complex plane by the field 
variable 0(a;^) as the coordinate x^ spans the entire space axis. Inequivalent vacuum 
configurations with winding numbers 0 and 1 respectively are depicted in Figs, la 
and Ic. In the contour of Fig. la the phase of 0 stays fixed at zero as x^ ranges 
between —cx) and +cx), whereas it goes once around the unit circle in Fig. Ic. Thus 
the corresponding vacuum configurations have winding number 0 and 1. The detailed 
variation of the phase is immaterial since it can always be changed locally by a gauge 
transformation. Thus, in Fig. la for example, as x varies from —cx) to +cx the 
field does not have to stay fixed, but could wander continuously on the unit circle 
provided the net change in phase is zero. However, the configuration of Fig. la 
cannot be continuously deformed to that of Fig. Ic without leaving the manifold of 
zero energy configurations. Therefore, in an evolution between neighboring vacua the 
field configuration must pass over an energy barrier, as illustrated in Fig. lb which 
singles out the configuration for which 0 vanishes at a point, rendering its phase 
there undefined. Figures Id-f add the additional perspective of spatial dependence 
for the field 0(a;^). Figures, la-c can be viewed as projections onto the complex plane 
orthogonal to the x^ axis of the curves in Figs. Id-e. 
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Figure 1: Example of two inequivalent vacuum configurations (a, c) and a field configuration at the 
top of the energy barrier separating them (b). Figures a, b and c trace the field (f) in the complex 
plane as the spatial coordinate spans the entire axis. A three dimensional perspective has been 
added in figures d, e and f to illustrate the detailed dependence of (j) on the spatial coordinate. 

The condition that 0 should take the same value at = ±cx) can be relaxed. 
Sometimes it is convenient to use the time independent gauge freedom to make 4>{oo) 
and 0(—cx)) differ (while keeping both hxed in time). Thus, the conhgurations of 
Figs, la-c can be gauge transformed into the conhgurations shown in Figs. 2a-c. In 
Fig. 2a the phase of 0 changes by —tt as goes from — cx) to +cx), whereas in Fig. 2c 
it rotates by tt. As in Fig. 1, the two vacuum conhgurations diher by a phase of 27r, 
i.e. by a unit change of winding number. In the intermediate conhguration (Fig. 2b) 
the held takes only imaginary values. In this gauge the conhguration which minimizes 
the energy on top of the barrier (i.e. the sphaleron conhguration) takes a very simple 
form 

0(a;^) = Ah[V^(a;^ — c)], = 0. (5) 

A possible parameterization for the entire evolution illustrated in Fig. 2 can be 
conveniently written 


0(a;^) 


1 — exp[iT — 2\/\{x^ — c) 

1 + exp[*r — 2V^(a;^ — c)] ’ 


Aq — 0 , Ai 


At^Tx 

7rch[2v/A(a:^ — c) 


( 6 ) 

(7) 
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Figure 2: A different gauge equivalent representation of the configurations illustrated in Fig. 1. 

As the reader can easily verify, for r = —7r/2 and r = 7r/2 the field 0 reduces to 
a number of unit modulus precisely spanning the contours of Fig. 2a and Fig. 2c 
respectively (as as ranges from —oo to +cxo). The corresponding values of Ai are 
chosen to make the gauge covariant derivative of 0 vanish, thus ensuring vacuum. 
We should point out, however, that Eqs. ||,|^ do not represent the solution of any 
special equations of motion (Euclidean or Minkowski). They are merely a compact 
parameterization of interpolating configurations, in terms of two variables c and r, 
which might be useful in studying sphaleron transitions based on the method of 
collective coordinates. 

If one couples chiral fermions to the gauge held in the 2-dimensional Abelian 
Higgs model the fermionic current has an anomaly which leads to fermion number 
violation in the topology changing processes described above. Thus this model would 
appear a very convenient system for a simplihed study of baryon number violation 
in high energy processes. However, as we will discuss in the next section, a crucial 
component of the computational investigation is the ability to numerically identify 
the normal mode amplitudes of the helds in the asymptotic regime. No matter how 
non-linear the system may be at any given point in its evolution, typically the energy 
will eventually disperse and bring the system to a regime where the helds undergo 
small oscillations about a vacuum conhguration. This dispersion is expected to occur 
in any held theoretical system, unless prevented by conservation laws such as those 
underlying soliton phenomena. Now, while the 2-dimensional Abelian Higgs model 
does not possess soliton solutions, we have observed computationally that the decay 
of the sphaleron in this system nevertheless gives origin to persistent, localized, large 
oscillations with an extremely small damping rate (this observation was also made by 
Arnold and McLerran in Ref. 0). These oscillations, illustrated in Fig. 3, make the 
system quite unwieldy for a computational investigation of baryon number violation 
based on semiclassical techniques. Thus we eventually turned our attention to the 
more realistic 4-dimensional SU (2) Higgs system, which constitutes the most relevant 
component of the full electroweak theory. Because of the larger dimensionality of 
space one would expect the energy to disperse much more readily in this system, an 
expectation borne out by results of Hellmund and Kripfganz||10||, who observed the 
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Figure 3: Sphaleron decay in the 2-dimensional Abelian Higgs model: evolution of the (j) field. The 
values of the phase of the complex field are coded by shades of gray, and the modulus of the field 
by the height of the surface. The sphaleron decays rather quickly, but leaves behind a quasi-stable 
oscillating remnant. 

onset of a linear regime following the sphaleron’s decay. 

The 3+1 dimensional SU{2) Higgs system is dehned in terms of a complex 
doublet <F(a;) and the gauge potential A^{x) with action 

S = Jdx^ - A(|$|2 - 1)2| , (8) 

where the indices run over 0 ■ ■ ■ 3 and where 

Ff,u = - dyA^ -%[A^, Ay ( 9 ) 

= (d^-zA^)^ ( 10 ) 

with Af^ = We use the standard metric = diag(l, —1, —1, —1), and have 

eliminated many inessential constants by a suitable choice of units. We focus on the 
spherically symmetric configurations of Ratra and Yaffe ||TT| , which reduce to an ef¬ 
fective 2-dimensional theory. This lower dimensional theory has the full topological 
structure of the 4-dimensional system, while having the virtue of being computation¬ 
ally manageable. 

The spherical Ansatz is given by expressing the gauge and Higgs fields in terms 
of six real functions Oq , Ui, a , /3, /i and z/ of r and t\ 


Ao(x,t) = -ao(r,t)cr-x 

Ai(x, t) = - [ai(r, t)a ■ x£* + - cr ■ xF) + 

2 r r 

$(x,t) = [/i(r,t)+ii/(r,t)o--x]^ , ( 11 ) 


where x is the unit three-vector in the radial direction and ^ is an arbitrary two- 
component complex unit vector. Note that conhgurations in the spherical Ansatz 
remain in the spherical Ansatz under gauge transformations of the form 

A^ ^A^ + iU^d^U /i = 0---3 (12) 

$ ^ f/$ , (13) 

where the gauge function is given by 

U = exp[*fl(r, f)cr ■ x/2] . (14) 


Inserting Eqs. 0 directly into Eq. one obtains an effective 2-dimensional theory 
with action 


S 


dvr 



- + (D'-xyD.x + 

“ Re(ix>^) 

-Ar^ . 


( 15 ) 


where the indices now run from 0 to 1 and in contrast to Ref. 
lowered with 77 ^,^ = diag(l, — 1 ), and where 


ir 


are raised and 


fdfj^tti/ Oi/ttfi ( 16 ) 

X = <y + ^P (17) 

0 = fi + tu (18) 

D^,x = {df,-iaf,)x (19) 

D^(t) = ^ «/ a )0 • ( 20 ) 

The reduced action, Eq. |I^, is invariant under the gauge transformation 

^ a^ + d^Vt ( 21 ) 

X ^ , (22) 

0 ^ (23) 
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which corresponds to the residual U{1) gauge invariance of Eqs. |^|T^. From Eqs. 


If® we see that the spherical Ansatz effectively yields a system very similar to 
the Abelian Higgs model considered above. In this reduced system the variables 
ao(r, t) and ai(r, t) play the role of the gauge held, whereas the variables x(r, f), which 
parameterizes the residual components of the 4-dimensional gauge held, and 0(r, f), 
which parameterizes the 4-dimensional Higgs held, both behave as 2-dimensional 
Higgs helds. Of course, the presence of metric factors (powers of r) in the action 
Eq. ^ is a reminder that we are really dealing with a 4-dimensional system. 

Regularity of the 4-dimensional held conhguration for r = 0 requires 


X{r = {)) = -I 

Im0(r = 0) = 0 . 


( 24 ) 


Although one could modify these boundary conditions by a singular gauge transfor¬ 
mation, this would introduce unnecessary complications. So we will impose Eqs. |2^, 
in addition to some other boundary conditions that also follow from the regularity of 
the 4-dimensional conhguration. It is also worth noting here that non-singular gauge 
transformations satisfy the condition that the gauge function f2(r, t) vanishes at r = 0 
(or is a multiple of 27r). 

We shall work in the Oq = 0 (or Aq = 0) gauge throughout. In the overlaying 
4-dimensional theory, if one compactihes 3-space to S'^ by identifying the points at 
inhnity, it is well known that the vacua correspond to the topologically inequivalent 
ways that S'^ can be mapped into SU{2) ~ These maps are characterized 


by the third homotopy group of SU{2), and a vacuum can be labeled by an integer 
called the homotopy index or winding number. The ehective 2-dimensional theory 
inherits a corresponding vacuum structure. From Eq. ^ it is apparent that the 
vacuum states are characterized by |x| = |0| = 1, with the additional constraint that 
(^-s well as Dix = Di(j) = 0). A convenient zero-winding vacuum is given 
by Xvac = — 4>vac = 1- Noutrivial vacua can be obtained from this vacuum via 
gauge transformations in which the gauge function H —>■ 27m (for non-zero integers 
n) as r —cx). Note that the compactihcation of 3-space ensures that at inhnity H 
is a multiple of 27r. Since H has been locked down to zero at r = 0, the winding of 
such a gauge transformed conhguration is just the integer n. For example, a typical 
winding-number one vacuum obtained from the previous trivial vacuum is given by 
Xvac = —^exp[^6'(r)], (jy^ac = exp[z6'(r)/2] and ai^ac = dr9{r), where 9{r) varies from 0 
to 271 as r changes from 0 to cxd. By taking advantage of the freedom of performing a 
time independent gauge transformation, however, one can also choose a gauge where 
x{r) and (j){r) tend to values diherent from —t and 1 as r —»• cx) (the condition 
IX*(t>^ —> — 1 must be preserved). Indeed, it will be convenient to choose one such 
gauge to parameterize the sphaleron. 

In making a topological transition between two inequivalent vacua, one must 
leave the manifold of vacuum conhgurations and pass over an energy barrier. Along 
such a trajectory there will be a conhguration of maximum energy. Of all these 
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Imij) 



(c) 


Ret]) 


Figure 4: The behavior of the (j) field for a typical topological transition. The X ii^is a behavior 
similar to the one in Fig. 2. 


maximal energy configurations, the sphaleron has the lowest energy and represents a 
saddle point along the energy ridge separating inequivalent vacua[§. The sphaleron 
can be expressed in the spherical Ansatz, and it is convenient to choose a gauge in 
which ttfj, = 0 and 


Xsph{r) = i[2f{r) - 1] (25) 

4>sph{r) = ih{r) , 

where / and h vary between 0 and 1 as r changes from 0 to oo and are chosen to 
minimize the energy functional. 

This choice of gauge for the sphaleron is slightly peculiar in the following sense. 
Finite energy conhgurations, like Eq. asymptote to pure gauge at spatial inhnity. 
Typically a gauge is chosen so that the appropriate gauge function is unity at spatial 
inhnity, and then space can be compactihed to the 3-sphere. But consider the spher¬ 
ical gauge function Eq. |T^ with f2(r) independent of time and varying from 0 to vr 
as r goes from 0 to oo. From Eqs. we see that the corresponding x, 0 —^ * for 

large r, which is the same as the sphaleron boundary conditions at spatial inhnity. 
Note, however, that U |r^oo= • x is direction dependent, so space cannot be com¬ 
pactihed. An arbitrary element of SU{2) can be parameterized by bol -b *cr ■ b where 
1 is the two by two unit matrix and b^ -|- = 1. Hence SU{2) ~ S^, and dehning 

the north and south poles by ±1, we see that ter ■ b with b^ = 1 parameterizes the 
equatorial sphere. Thus, the sphere at inhnity is mapped onto the equatorial sphere 
of SU{2). In this gauge, a topology changing transition proceeding over the sphaleron 
corresponds to a transition where the helds wind over the lower hemisphere of SU{2) 
before the transition and over the upper hemisphere after the transition, with a net 
change in winding number still equal to one. In this gauge, the behavior of the y 
held in a topological transition will be very similar to the behavior of the Higgs held 
in the 2-dimensional model, already illustrated in Fig. 2. The behavior of the 0 held 
is illustrated in Fig. 4. We could of course work in a compactihed gauge where a 
topological transition would occur between a held with no winding and a held with 
unit winding, as in Fig. 1, but the sphaleron would look more complicated. The 
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Figure 5: Behavior of the X in a topology changing transition. The various shades of gray 
code the phase of the complex field. The field starts as an excitation about the trivial vacuum, 
passes over the sphaleron and then emerges as an excitation about the vacuum of unit winding. 
Note the persistent strip of 27r phase change after the wave bounces off the origin. 


advantage of Eq. ^ from a computational perspective is that perturbations about 
the sphaleron can be more easily parameterized. 

In the sphaleron configuration the y held has a zero at some non-zero value of 
r, whereas the (j) held has a zero for r = 0 corresponding to the vanishing at the origin 
of the actual 4-dimensional Higgs held. The zero of y is reminiscent of the zero which 
characterizes the sphaleron of the Abelian Higgs model. However, as shown in Ref. 

, it is the zero of the true Higgs held (i.e. the zero of 0) which carries a deeper 
signihcance and should be associated with the actual occurrence of the topological 
transition. Nevertheless, since the phase changes of x are more dramatic than those of 
0, for purposes of illustration it is often more convenient to plot y. Figure 5 illustrates 
a typical topological transition. The conhguration starts as a small excitation about 
the trivial vacuum dehned above, passes over the sphaleron and then emerges as a 
conhguration that undergoes a 27i phase change. 

3 Numerical study of classically allowed processes 

For energy larger than the sphaleron energy Egph, i-e. for e > egph = g'^Egph, 
classical evolution which changes the topology of the helds becomes possible. So¬ 
lutions to the classical equations of motion provide the dominant contributions in a 
weak coupling expansion of the path integrals describing transition processes between 
coherent states. The existence of topology-changing solutions indicates that the cor- 



12 




responding processes, which because of the change of topology violate baryon number, 
are not suppressed. However it would be premature to conclude that baryon number 
violation can occur with non-negligible amplitude in high energy collisions. Indeed, 
the number of particles in the initial state of such processes is small and, in terms of 
the rescaled helds used in the description of the classical evolution, this converts into 
a rescaled particle number v = which tends to zero as g —>■ 0. Thus, in order 
to establish evidence for baryon number non-conservation in high energy collisions, 
one must show that topology changing classical evolution can occur with arbitrarily 
small u in the initial state. 

The primary impediment for rapid baryon number violation seems to be a phase 
space mismatch between initial states of low multiplicity and hnal states of many 
particles. The authors of Ref. |T^ look at simplihed models and observe that, clas¬ 


sically, it is difficult to transfer energy from a small number of hard modes to a large 
number of soft modes. However, Ref. [|T^ hnds that for pure Yang-Mills theory in 


2 -dimensions the momenta can be dramatically redistributed, but unfortunately the 
initial particle number seems to be rather large in their domain of applicability. It is 
the purpose of our investigation to shed light on the situation in 4-dimensions when 
the Higgs held is added. 

Given any classical evolution, because of the dispersion of the energy, the helds 
will asymptotically approach vacuum values. Thus, for times t < —Ti and t > Tf and 
sufficiently large T*, Tf, the equations of motion will reduce to linearized equations 
describing the small oscillations of the system about a vacuum conhguration. In 
this linear regime, the evolution of the helds will be given by a superposition of 
independent harmonic oscillators (the normal modes). In terms of the frequencies a;„ 
and amplitudes a„ of these oscillators the (rescaled) energy and particle number are 
given by 

^ (26) 

(27) 




and 


U — ^(0 Cln<ln- 


Thus for any classical evolution the energy e and the particle numbers z/j and Uf 
of the asymptotic initial and hnal states are well dehned. (The energy is of course 
conserved and well dehned even in the non-linear regime.) In addition, because of 
the fact that the helds approach vacuum values for t —>■ ±oo, the winding numbers of 
initial and hnal conhguration are also well dehned. Because of the sphaleron barrier, 
the energy e of all the classical solutions with a net change of winding number is 
bounded below by the sphaleron energy egph- The problem one would like to solve, 
then, is whether the initial particle number z/j of these solutions can be arbitrarily 
small, or more generally, one would like to map the region spanned by all possible 
values of e and Ui for the topology changing classical evolution. The highly non-linear 
nature of the equations of motion makes an analytic solution unlikely, even if one is 
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willing to make the crudest approximations. The problem is however amenable to 
solution by computational techniques. In this section we will illustrate the strategy 
we have followed to formulate it on the computer and the progress we have been able 
to make towards its solution. 

The fundamental computational ingredient consists in the implementation of 
a numerical solution of the equations of motion. We start from the Hamiltonian 
formulation of the equations of motion for the continuum system in the oq = 0 gauge. 
In such a formulation the variables 


ai(r), x(r), 0(r) 


form a set of canonical coordinates, conjugate to the momenta 

E{r) = r'^dotti, 

^xi^) = 9ox, 

71-0(r) = r^(9o0 . 


The evolution of these variables is generated by the Hamiltonian 


H = 



^ + kxP + + \DrX\^ + r^\Dril>f 

+Ar^(0p-l)^ . 


(28) 


(29) 


(30) 


Gauss’ law 

drE = - X*T^x) + - ^*^ 4 ) ( 31 ) 

expresses the residual invariance of the system under time independent local gauge 
transformations and is imposed as a condition on the initial state. It is then auto¬ 
matically conserved by the equations of motion. 

To solve the equations of motion numerically the system must be discretized. It 
is convenient to use the formalism of lattice gauge theories. The r-axis is subdivided 
into N equal subintervals of length Ar with hnite length L = N Ar. Thus, the lattice 
sites have spatial coordinates r* = iAr with i = 0 ■ ■ ■ N, and the midpoints between 
lattice sites have coordinates rj+ 1/2 = (f -|- 1/2) Ar with i = 0 ■ ■ ■ N — 1. The fields y, 
0, TT^ and 7r0 are then represented by discrete variables dehned over the end points of 
the intervals, Xi = (t^i = 0(a), etc.; whereas the helds ai and E are dehned over 

the intervals themselves by an = ai = ai(ri+i/ 2 ) and Ei = E{ri+i/ 2 ) (for notational 
simplicity we have dropped the spatial subscript on the discretized gauge held). The 
covariant derivatives are then replaced by covariant hnite diherences, e.g. 


DrX 


exp[-mj Ar] Xi+i - Xi 
Ar 


i = 0---A- 1 , 


(32) 
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and like the gauge fields they are to be thought of as living on the links between lattice 
sites. The rest of the discretization is straightforward. One obtains a discretized 
Hamiltonian expressed in terms of a finite set of variables, which still possesses 
exact local gauge invariance under the transformations of Eqs. provided that 

Oi —> Oi + dr^l is replaced with the finite difference formula 


n 


CLi 






Ar 


i = 0...N-l , 


(33) 


where Oj = 0(rj). From Hj:, one can easily obtain the canonical evolution equations 
for the discretized variables. Gauss’ law, which now takes a discretized form, must 
be imposed on the initial state and is then preserved (exactly) by the time evolution 
because of the gauge invariance of the discretized system. In practice we have used 
values of N equal to 256, 512 and 1024 and values of Ar equal to 0.2, 0.1 and 0.05 
in respectively study the properties of the system (with A = 0.1). We found these 
parameters to be adequate for obtaining, on the one hand, a reasonable approximation 
to the continuum system and, on the other, a cut-off on r sufficiently large to allow 
for an effective linearization of the equations of motion before the waves hit the 
boundary. The restriction to uniform spacing of the subintervals on the r-axis is not 
fundamental and we have also implemented a discretization where Ar increases as 
one moves out on the r-axis. In this manner one can effectively make the system 
larger and delay the effects of the impact of the waves with the boundary without 
worsening the spatial resolution near r = 0, where most of the non-linear dynamics 
takes place. We have found however that the advantages one gains hardly warrant 
the additional complications introduced by the non-uniform spacing. 

For the numerical integration of the time evolution we have used the leap-frog 
algorithm. Since this algorithm, or the equivalent velocity Verlet algorithm, consti¬ 
tutes one of the fundamental techniques for the integration of ordinary differential 
equations of the Hamiltonian type and as such is textbook material, we will not dis¬ 
cuss it in any great detail. Essentially, given conjugate canonical variables g* and Pi 
which obey equations 


dqi 


dpi 

dt 


= fi{q), 


(34) 


one evolves the values of q and p from some initial t to t + At as follows. In a first 
step Pi is evolved to the mid-point of the time interval by 


Pi 

qi 


Pi= Pi + fM)^^ 


Qi = 


qi 


(35) 


(Although qi is left unchanged, it is convenient to consider the step formally as a 
transformation of the entire set of canonical variables.) In a second step one evolves 
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the coordinates from their initial value g* = q[ to their value at the end of the interval 


P’i 


Pi = 
Qi = 


P'i, 


q'i +■ 


(36) 


Finally the momenta are evolved from their value at the midpoint to the hnal value 


p" ^ Pf = P" + /i(9")^. 

Qi g” = g” ■ 


(37) 


One can easily verify that these equations reproduce the correct continuum 
evolution from f to f + At up to errors of order (Af)^. Moreover, the algorithm has 
the very nice property that all three steps above constitute canonical transformations 
and that it is reversible (in the sense that starting from g'", —p'-', up to numerical 
errors one would end up exactly with g*, —pi). Another very nice feature of the 
algorithm is that, although the evolution of the variables is affected by errors of order 
(At)^, the energy of a harmonic oscillator, and therefore also of any system which 
can be decomposed into a linear superposition of harmonic oscillators, is conserved 
exactly (always up to numerical errors). In a sequence of several iterations of the 
algorithm after the momenta have been been evolved by the initial At/2, the hrst 
and third step can be combined into a single step whereby the momenta are evolved 
from the midpoint of one interval to the midpoint of the next one, “hopping over” 
the coordinates, which are evolved from endpoint to endpoint. This motivates the 
name assigned to the algorithm. 

With a good grasp on the numerical solution of the equations of motion, we can 
now turn to the second crucial component of the computation, namely the identih- 
cation of the particle number in the initial state. One could easily parameterize an 
initial conhguration of the system consisting of incoming waves in the linear regime. 
However, it would be extremely difficult to adjust the parameters so as to insure that 
a change of winding number occurs in the course of the subsequent evolution. For 
this reason it is much better to parameterize the initial conhguration of the system 
at the moment when a change of topology occurs. Thus our strategy consists in 
implementing a time-reversed solution of the equations of motion, where the initial 
conhguration is the conhguration of the system at the moment when it passes over 
the sphaleron barrier and the asymptotic conhguration for large t will be interpreted 
as a time-reversed incoming state. 

Topology changing transitions within the spherical ansatz are characterized by 
the vanishing of 0 at r = 0 and the the vanishing of y at nonzero r. For a sequence 
of conhgurations that pass directly through the sphaleron these two zeros occur at 
the same time. However, this is not the most general case and the zeros of 0 and y 
need not occur simultaneously |^. So we have parameterized initial conhgurations 
in terms of coefficients Cn of some suitable expansion of the helds and their conjugate 
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momenta, constrained only by the fact that the held y must have a zero at some 
hnite r, as in the sphaleron conhguration. Furthermore, we can use gauge invariance 
to make this held pure imaginary. The held 0 is only restricted to obey the boundary 
conditions and does not necessarily vanish at the 
some point in its evolution). 0 

In practice, computational considerations will limit the number of parameters 
we will be able to use. We will then evolve the system until the dispersion of the 
energy brings it to the linear domain. At this point, we can calculate the ampli¬ 
tudes of the normal modes of oscillation and the particle number of the system in its 
asymptotic state. It is clear that every set of coefficients c„ will determine one dehnite 
value for the energy e and particle number v of the asymptotic state. In calculating 
the particle number, one should use only the lower lying modes since higher modes 
probe wave lengths of order the lattice spacing. For our lattice parameters, we found 
that considering the hrst 50 to 100 modes is reasonable. We should point out that 
an arbitrary initial conhguration is not necessarily guaranteed to change topology. 
However, by evolving the conhguration both forward and backward in time we can 
easily verify whether topology changes, and initial conhgurations that do not change 
the topology can be rejected. By reversing the time evolution, then, we will have de- 
hned an initial asymptotic conhguration with energy e and particle number i/, which 
in the course of its evolution undergoes a change of topology. By varying the values 
of the parameters we will be able to study the properties of such held evolution 
and, in particular, explore the domain of permissible values for e and v. It should 
be obvious at this point that the determination of the normal mode amplitudes is 
another crucial ingredient of the computation. 

The normal modes of oscillation can be calculated starting from the linearized 
equations of motion. These equations can be obtained from an expansion of the 
Hamiltonian up to second order in the deviation of helds from a vacuum conhguration. 
For conciseness of presentation we will not reproduce here their explicit form, but will 
limit ourselves to a discussion of the general properties of the normal modes. The 
normal modes can be obtained by assuming an oscillatory evolution for the helds of 
the type 


origin (although it will vanish at 


sm[uJnt] 

cos[uJnt] . (38) 

In this equation we have made reference only to one pair of conjugate helds, but the 
equation should be complemented with a similar assumption for the pairs Oi and E, 
6(j) and tt,^, which are a priori all coupled together with 6x and in the equations 
of evolution. We have used the symbol 6 to denote the deviations of the helds y and 
0 from their vacuum values x = —* and 0=1. All of the other helds vanish in the 

“Of course we could equally well arrange (j) to be pure imaginary and to vanish at the origin, with 
no restriction y on other than its boundary conditions. 


17 



trivial vacuum configuration. For the discretized system, r should be replaced by an 
index i. 


Substituting the Ansatzoi Eqs. ^into the equations of motion one obtains a set 
of eigenvalue equations. Their solution determines the possible values of ojn as well 
as the corresponding eigenfunctions (or more properly eigenvectors in the discretized 
case) 'K^\r) etc. We have determined eigenvalues and eigenfunctions both 

numerically, for the discretized system, and analytically, for the continuum system. 
One hnds that the normal modes of oscillations are naturally grouped together into 
four sets of modes: 

i) a set of modes where only the imaginary part of the helds 5x and tt^ is 
non-vanishing. These correspond to an oscillation of the modulus of the y held . 

ii) a set of modes where only the real part of the helds 50 and tt^ is non-vanishing. 
These correspond to an oscillation of the modulus of the 0 held . 

hi) and iv) two sets of modes where the real part of 5%, vr^, the imaginary 
part of 50, as well as oi and E are coupled together and non-vanishing. These 
correspond to oscillations of the phases of y and 0 (in coherence or opposition of 
phase), accompanied corresponding oscillations of the gauge helds. 

Given the expressions for the eigenfunctions it is possible to extract the ampli¬ 
tudes of the normal modes of oscillation an{t) by taking suitable convolutions of the 
eigenfunctions with the helds and momenta. This procedure exploits various prop¬ 
erties of orthogonality which the eigenfunctions satisfy. One subtle point, however, 
involves the need to hx the gauge. The normal modes are obtained on the basis of an 
expansion into small oscillations around the trivial vacuum conhguration with con¬ 
stant helds (one could of course expand around any hxed vacuum conhguration, but 
the formulae are much simpler if one expands around constant helds). However there 
is no guarantee that the evolution of the helds will lead to an asymptotic regime of 
small huctuations precisely around such a vacuum conhguration. Indeed, in general 
this will not happen and the actual conhguration will typically diher from the one 
used to derive the normal modes by a large gauge transformation. The remedy is 
easy enough. One can perform a gauge transformation to a hxed gauge which dihers 
at most by small huctuations from the one where the expansion has been performed. 
In our computations we have used the unitary gauge dehned by Argx(r) = —7r/2, 
transforming the helds to this gauge at some dehnite time Iq, at which point the 
appropriate amplitudes may be extracted. At subsequent times the helds might no 
longer be small perturbations about the trivial vacuum, so every time the amplitudes 
are calculated we hrst enter the above gauge. An alternative approach, which we 
have also implemented, consists in deriving the linearized equations of motion for 
a complete set of gauge invariant quantities |^(we have used |x| — 1, |0| — 1, the 
electric held E, the diherence of phases Argy — 2Arg0 and the time derivatives of 
all these quantities). Following a procedure similar to the one outlined above (cf. 
Eq. |3^) one can derive the eigenvalues and eigenfunctions for the small oscillations 
of these quantities (the eigenvalues are identical, of course, to those obtained using 
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Figure 6: Sphaleron decay in the four dimensional SU(2) Higgs model: evolution of the X field. 
The values of the phase of the complex field are coded by different shades of gray. When the evolution 
reaches the linear regime, a gauge transformation, indicated by the sudden change of shading, is 
performed to extract the normal mode amplitudes. 


gauge variant quantities) and extract the amplitudes through suitable convolutions 
with the evolving helds. We have verihed that the two procedures produce identical 
results. 

We can now turn our attention to the hgures that illustrate our results for a 
system with coupling constant A = 0.1. Figure 6 illustrates the behavior of the held x 
following the decay of the sphaleron after a slight initial perturbation. We have found 
it very convenient and informative to use color to code the phase of the complex helds. 
Unfortunately the illustrations in these pages cannot be reproduced in color and we 
have tried to render the variation of the phase with a gray scale. From Fig. 6 it is clear 
that the energy, which is concentrated in the neighborhood of r = 0 in the sphaleron 
conhguration, disperses and gives rise to a pattern of outgoing waves. The sudden 
variation of tonality at some point in the time evolution indicates the change of phase 
induced by the gauge transformation to the unitary gauge. In Fig. 7 we display the 
behavior of the particle number in the four normal modes of oscillations as function of 
time. It is apparent that, after an initial transition period where the particle nnmber is 
not well dehned, the qnantities settle to valnes which are reasonably constant in time. 
We take this as evidence that the system has reached an asymptotic regime where one 
can meaniningfnlly dehne a conserved particle nnmber. Finally, the evolntion of the 
system can be time reversed, as we have discnssed above, and the (time-reversed) hnal 
conhgnration can be considered as an asymptotic initial conhguration with dehnite 
energy and particle nnmber. The previously shown Fig. 5 illnstrates the evolntion 
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Figure 7: Sphaleron decay in the four dimensional SU(2) Higgs model: behavior of the particle 
number in the four normal modes of oscillation of the linearized system as function of time. The 
sum in Eq. 27 extends over the first 50 modes for each of the four mode types. 


obtained taking one such asymptotic initial configuration. The incoming waves are 
seen to merge in the neighborhood of the origin, where a change of topology takes 
place. The fact that the winding number of the field configuration has changed is 
indicated by the strip of rapidly varying tonality which persists in the neighborhood 
of the origin and codes the variation of the phase of y. With color, this strip would 
appear as a vivid rainbow, left over as a marker of the change of topology of the 
evolving fields. 

The configuration which sits on top of the sphaleron barrier can be parameter¬ 
ized by expanding into suitable components a complete set of independent fields. We 
have chosen these fields to be the perturbations Ay and A(j) (not necessarily small) 
of the y, (f) fields of the sphaleron, the field ai and the momenta tt^, We avoid a 
redundant gauge degree of freedom by taking Ay purely imaginary (like the y field 
of the sphaleron itself). The final field needed to specify the initial configuration, 
i.e. the electric field E, can then be derived from Gauss law. We have parameterized 
these fields in terms of Bessel functions, chosen in such a way to respect appropriate 
boundary conditions at r = 0 and r = L. The coefficients c„ of the expansion can 
now be varied and, in this way, one can explore the region in the e-u plane spanned 
by all of the topology changing classical solutions. 

Of course, the space of topology changing configurations is infinite and a random 
exploration of such space would not lead to very useful results. We must keep in mind 
that the interesting question is whether there is a lower bound in u or, more generally, 
what is the lower boundary of the region spanned by all topology changing solutions 
in the e-u plane. This question can be investigated by methods of stochastic sampling. 
One can perform random small steps in the space of all configurations by varying the 
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Figure 8: Monte Carlo results with lattice parameters of Ar = 0.2, N = 256 and L=51.2 and 
Higgs coupling A = 0.1. 

parameters c„ stochastically. After each change the new conhguration is evolved until 
one can extract the particle number in the linear regime. Then the variations of energy 
and particle number Ae and Au induced by the change Ac^ become well dehned. A 
standard Metropolis Monte Carlo sampling technique consists in accepting or rejecting 
the change according to the value of the quantity AF = pAe + /iAz/, where P and 
fi are parameters that weight what region of the e-u plane is explored. To be more 
precise, the change is accepted with conditional probability P = Min[l,exp(—AF)], 
which has the effect of producing configurations distributed according to a measure 
proportional to exp(—/3e — /iz/). The parameters P and /i play the role of inverse 
temperature and chemical potential. By choosing these parameters appropriately 
one can drive the sampling towards low values of energy and particle number and 
thus explore the interesting region in the e-v plane. We have begun implementing 
this procedure and Fig. 8 illustrates our hrst results. It is interesting to note that the 
decay of a (slightly perturbed) sphaleron gives rise to a particle number v k, 1.9(47r). 
For g = 0.6 this corresponds to N k, QQ physical particles. From Fig. 8 one can see 
that our sampling procedure has produced conhgurations with comparable energy 
and much smaller particle number. Of course the ultraviolet cutoff induced by the 
lattice discretization puts a lower limit on the ratio z//e, which occurs when only the 
highest mode is excited. We have used lattice parameters Ar = 0.2, N = 256 with 
L = 51.2 and have considered only the hrst 50 modes for each of the four types of 
normal modes. Hence the minimum value of z//e is of order l/oJmax ~ L/umaxT^ ~ 0.3. 
Given our lattice resolution, we have saturated the lowest bound in particle number 
that we are sensitive to. This may be an indication that there is no lower bound on 
z/, but our calculation is still at a very preliminary stage and much more work will be 
needed to establish reliable results. 
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Figure 9: Complex time contour. 

4 High energy baryon number violating processes 
below the sphaleron barrier 

In the work of Ref. 0 Rubakov, Son and Tinyakov relate the probability of 
a topology changing process from a coherent state to the action of a complexihed 
classical evolution along a special contour in the complex time plane. As shown in 
Fig. 9, this contour consists of a semiaxis parallel to the real axis a = {—oo + tT —>■ 
tT), followed by a segment (3 = {iT —0) of length T along the imaginary axis, 
followed by the real positive semiaxis 7 = (0 —>■ +cxd). The word “complexihed 
classical evolution” refers to the fact that the equations of motion must be analytically 
continued to complex helds. This has to be handled with some care for the case in 
which helds, such as x(r) and 0(r), are already complex to begin with. In this case 
x(r) and its complex conjugate, which we will denote here by must be considered 
as formally independent variables, which can be analytically continued separately, and 
do not necessarily satisfy the “reality condition” x(r) = x(r)*. The same applies to 
0 (r) and ^(r). 

The boundary conditions are of special importance. The solution must be real, 
in the sense specihed above, along 7 , where the asymptotic helds for t +cx) repre¬ 
sent the hnal state. The condition on the particle number in the initial state translates 
instead into the requirement that for very early times in the evolution (i.e. asymptot¬ 
ically ioi t ^ —00 + iT along a) the helds must reduce to a superposition of normal 
modes of oscillation with amplitudes satisfying the equation 

a{—k) = e^a*{k) (39) 

where k is an index characterizing the radial momentum of the waves, and 9 plays 


22 





the role of a chemical potential conjugate to the particle number in the initial state. 
It is clear from Eq. ^ that for 6 * 7 ^ 0 the helds cannot be real along a. 

Computationally, it is convenient to think again of a time reversed evolution by 
which, starting from an initial conhguration at t = 0 the fields undergo a Euclidean 
evolution along the imaginary time axis to t = iT (i.e. following the oriented segment 
—f3) and then a Minkowski evolution (along —a) from tT to —00 + tT. For t = 0 
for an a:-axis discretized with N sites we have as free variables N real values of the 
field and N pure imaginary values of the conjugate momentum per each independent 
canonically conjugate held-momentum pair. The conditions that the normal mode 
amplitudes must satisfy Eq. ^ amount to N complex equations, i.e. 2N real equations 
again for each canonically conjugate held-momentum pair. Thus, in principle, one 
could evolve the helds from an initial Ansatzat t = 0 and adjust the initial variables so 
that Eq. ^ is satished. In practice, since the evolution equations along the imaginary 
time axis are elliptic, one cannot perform a forward integration. Rather, one must 
resort to some relaxation procedure or other global algorithm, by which one solves 
the evolution equations as a set of simultaneous non-linear equations for all points 
of a space-time grid. The situation is further complicated by the fact that, with 
complexihed helds, one cannot just minimize a Euclidean action integral. We have 
developed a “second order” formulation, by which we minimize a constraint functional 
obtained from the modulus squared of the functions that must vanish at all grid points 
(an earlier part of this study was done in collaboration with Timothy Vaughan). In 
this case also we have used the formalism of lattice gauge theory to obtain a gauge 
invariant discretization. We tested our procedure in the context of the 2D Abelian 
Higgs model (one space, one time dimensions), where we found that it did reproduce 
the expected Euclidean solutions, including solutions with multiple bounces of the 
helds between two diherent topological sectors. 

Another crucial component of the calculation consists in solving the evolution 
equations along the semiaxis —a, from tT to —oo + tT and extracting the normal 
mode amplitudes a{k),a{k). Here we feel that we have formalism and algorithm 
already in place, although we need to check that the integration remains stable with 
complexihed helds. 

In conclusion, the study of topology changing processes below the sphaleron 
barrier presents some additional challenges, the most notable being the need to inte¬ 
grate elliptic, rather than hyperbolic, equations over part of complex time path. The 
number of degrees of freedom and the complexity of the calculation, however, are 
not substantially diherent from those which characterize the numerical investigations 
of classically allowed processes, where, as we have shown above, the power of the 
computational tools are well adequate to produce interesting and accurate results. 
This warrants the expectation that even the classically forbidden processes will be 
amenable to a successful computational analysis. 
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